Valley-Polarized Metals and Quantum Anomalous Hall Effect in Silicene 
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Silicene is a monolayer of silicon atoms forming a two-dimensional honeycomb lattice, which shares almost 
every remarkable property with graphene. The low energy structure of silicene is described by Dirac electrons 
with relatively large spin-orbit interactions due to its buckled structure. The key observation is that the band 
structure is controllable by applying electric field to silicene. We explore the phase diagram of silicene together 
with exchange field M and by applying electric field Ez. There appear quantum anomalous Hall (QAH) in- 
sulator, valley polarized metal (VPM), marginal valley polarized metal (M-VPM), quantum spin Hall (QSH) 
insulator and band insulator (BI). They are characterized by the Chem numbers and/or by the edge modes of a 
nanoribbon. It is intriguing that electrons have been moved from a conduction band at the K point to a valence 
band at the K' point for Ez > G in the VPM. We find in the QAH phase that almost flat gapless edge modes 
emerge and that spins form a momentum-space skyrmion to yield the Chern number. It is remarkable that a 
topological quantum phase transition can be induced simply by changing electric field in a single silicene sheet. 



Silicene, a monolayer of silicon atoms forming a two- 
dimensional honeycomb lattice, has been synthesized|tl|-|31] 
and attracts much attention |J4|48|] recently. Almost every strik- 
ing property of graphene could be transferred to this innova- 
tive material. It has additionally a salient feature, that is a 
buckled structure||4l|5|] owing to a large ionic radius of silicon. 
Silicene has a relatively large spin-orbit (SO) gap of 1.55meV, 
which provides a mass to Dirac electrons. Furthermore, we 
may control experimentally the massf?] by applying the elec- 
tric field Ez- Silicene undergoes a topological phase transition 
from a quantum spin Hall (QSH) state to a band insulator (BI) 
as I^E^I increases|7]. A QSH state is characterized by a full 
insulating gap in the bulk and helical gapless edgesll9l 4l2ll . 

There exits another state of matter in graphenel] 
is a quantum anomalous Hall (QAH) state 1 1 81 



llMJ, that 
1911 . charac- 
terized by a full insulating gap in the bulk and chiral gap- 
less edges. Unlike the quantum Hall effect, which arises from 
Landau-level quantization in a strong magnetic field, the QAH 
effect is induced by internal magnetization and SO coupling. 
In this paper we analyze the band structure of silicene to- 
gether with exchange field M and by applying electric field 
Ez to silicene. We explore the phase diagram in the Ez-M 
plane. Silicene has a rich varieties of phases because the elec- 
tric field Ez and the exchange field M have different effects 
on the conduction and valence bands characterized by the spin 
and valley indices. There are insulator phases, which are the 
QSH, QAH and BI phases. There emerges a new type of 
metal phase, the valley-polarized metal (VPM) phase, where 
electrons have been moved from a conduction band at the K 
point to a valence band at the K' point for Ez > 0. Such a 
phase is utterly unknown in literature as far as we are aware of. 
There are also metallic states on phase boundaries, which are 
metal (M), marginal- VPM (M-VPM) and spin VPM (SVPM) 
states. All these phases and states are characterized by the 
Chern numbers and/or by the edge modes of a nanoribbon. It 
is possible to materialize any one of them by controlling Ez at 
an appropriate value of M. Furthermore, as we have pointed 
out elsewherelLZl], by applying an inhomogeneous field Ez, it 
is possible to materialize some of these topological phases to- 
gether with states on the phase boundaries simultaneously in 



a single silicene sheet. 

Silicene consists of a honeycomb lattice of silicon atoms 
with two sublattices made of A sites and B sites. The states 
near the Fermi energy are tt orbitals residing near the K and 
K' points at opposite corners of the hexagonal Brillouin zone. 
We refer to the K or K' point also as the K^ point with the val- 
ley index 77 = ±1. We take a silicene sheet on the a;y-plane, 
and apply the electric field Ez perpendicular to the plane. Due 
to the buckled structure the two sublattice planes are separated 
by a distance, which we denote by 2£ with £ = 0.23A. It gen- 
erates a staggered sublattice potential a 2£Ez between silicon 
atoms at A sites and B sites. 

The silicene system is described by the four-band second- 
nearest-neighbor tight binding model, 
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where cj^ creates an electron with spin polarization a at site i, 
and {i,j) / ((«,.?)) i"un over all the nearest/next-nearest neigh- 
bor hopping sites. We explain each term, (i) The first term rep- 
resents the usual nearest-neighbor hopping with the transfer 
energy t ~ 1. 6eV. (ii) The second term represents the effective 



SO coupling with Aso — 3.9meV, where cr 
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the Pauli matrix of spin, with Vij = +1 if the next-nearest- 
neighboring hopping is anticlockwise and i/ij = —1 if it 
is clockwise with respect to the positive z axis, (iii) The 
third term represents the first Rashba SO coupling associ- 
ated with the nearest neighbor hopping, which is induced by 
external electric field[14, 20]. It satisfies Ari(0) = and 
becomes of the order of lO/ieV at the critical electiic field 
Ec = Aso/^ = 17meVA^^. (iv) The forth term represents the 
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FIG. 1: (Color online) Phase diagram in the E^-M plane. Heavy 
lines represent phase boundaries, where the system becomes metal- 
lic. Chem and spin-Chem numbers (C,Ca) are well defined and 
given in insulator phases. Dotted lines represent the points where the 
band gap closes, which are within the VPM phase. A circle shows a 
point where the energy spectrum is calculated and shown in Fig|2] 



second Rashba SO coupling with Ar2 = 0.7meV associated 
with the next-nearest neighbor hopping term, where /i^ = ±1 
for the A (B) site, and d.^ = dy/ |dy | with the vector dij 
connecting two sites i and j in the same sublattice. (v) The 
fifth term is the staggered sublattice potential term, (vi) The 
sixth term represents the exchange magnetization: Exchange 
field M may arise due to proximity coupling to a ferromag- 
net such as depositing Fe atoms to the silicene surface or de- 
positing silicene to a ferromagnetic insulating substrate, as has 
been argued for graphene| 13-15]. The Hamiltonian dl) can 
also be used to describe germanene, which is a honeycomb 
structure of germanium] 5, 6], where various parameters are 
t = 1.3eV, Aso = 43meV, Ar2 = lO.TmeV and i = 0.33A. 

In this paper we derive the topological phase diagram in 
the Ez-M plane and make its physical interpretation. The 
topological quantum numbers are the Chem number C and 
the Z2 index. If the spin Sz is a good quantum number, the 
Z2 index is identical to the spin-Chem number Cg. They are 
defined when the state is gapped and when the Fermi level 
is taken within the gap, and given by C = C+ + C_ and 
Cs — |(C'+ — C_), where C± is the summation of the Berry 
curvature in momentum space over all occupied states of elec- 
trons with Sz ~ ±1. They are well defined even if the spin is 
not a good quantum numberr 151 12111 . In the present model the 
spin is not a good quantum number because of spin mixing 
due to the Rashba couplings Ari and Ar2, and the resulting 
angular momentum eigenstates are indexed by the spin chiral- 
ity s = ±1. We can calculate these numbers at each point in 
the Ez-M plane by using the standard formulas lll3l - ll5ll . 




We present our result on the phase diagram in FigJT] We 
show later how to deiive the phase boundaries based on the 
low-energy Dirac theory. We have also calculated the band 
structure of a silicene nanoribbon with zigzag edges, which 
we give in Fig |2] for typical points in the phase diagram. The 



FIG. 2: (Color online) The band structure of a silicene nanoribbon at 
marked points in the phase diagram (Fig[T}. The vertical axis is the 
energy in uit of t, and the horizontal axis is the momentum. We can 
clearly see the Dirac cones representing the energy spectrum of the 
bulk, as explained in Fig[3] Lines connecting the two Dirac cones 
are edge modes. The spin s^ is practically a good quantum number, 
which we have assigned to the Dirac cones. An enlaraged figure of 
the QAH state is given in Fig[3] 



topological numbers are {C,Cs) — (0,0) in the BI phase, 
(0, 1) in the QSH phase, (2, 0) in the QAH phase with A/ > 
and (-2, 0) in the QAH phase with M < 0. In all these states 
the band gap is open, where the Fermi level is present, and 
they are insulators. 

We first discuss the system at £^2 = and compare our re- 
sults with those previously obtained in graphene O 131 - 1 1511 . The 
main difference is the appearence of almost flat edge modes 
in our system (Fig|3]i. This occurs because the Rashba in- 
teractions are different between these two systems. We have 
Ari — for Ez = and Ar2 = at the K and K' points in 
silicene, but Ari ^ and Ar2 = in graphene. Nevertheless, 
the difference is only quantitative. As far as the topological 
properties are concerned, there exists no difference. Indeed, 
in these two systems, the Chern number is identical in each 
corresponding phase together with quantized Hall conductiv- 
ity, and the edge states support the edge current. However, the 
group velocity of the edge modes is extremely small due to 
the almost flat gapless modes in silicene. 

Our most important result is the VPM phase, which appears 
in such regions that EzM ^ and occupies a major part of 
the phase diagram. A part of the conduction (valence) band is 
above (below) the Fermi level at the K (K') point for Ez > 0, 
as is observed in Fig|2VPM). Hence, electrons are moved 
from the K valley to the K' valley, as implies the valley po- 
larization. The phase is characterized by the property that it 
is a metallic state though gaps are open both at the K and K' 




FIG. 3: (Color online) The energy spectrum of a QAH state. Gray 
curves are for a nanoribbon, which are identical to Fig|2]QAH). 
Heavy green curves represent the energy spectrum of the bulk, calcu- 
lated independently. A gap opens in the bulk spectrum, where almost 
flat gapless modes appear at the edges of a nanoribbon. A red (blue) 
arrow indicates the spin direction away from the Fermi level. 



points. We note that the Chern and spin-Chem numbers are 
ill-defined in the VPM phase, since the Fermi level does not 
lie inside the band gaps at the K and K' points simultaneously. 

There exist M-VPM states on phase boundaries indicated 
by heavy lines in the phase diagram, where the conduction 
and valence bands touch the Fermi surface at the K and K' 
points, respectively, for Ez > 0. On the other hand, in SVPM 
states the conduction and valence bands touch the Fermi sur- 
face both at the K and K' points. We expect topological quan- 
tum critical phenomena in these states. 

In order to explore the physics underlying the phase dia- 
gram, we analyze the low-energy effective Hamiltonian de- 
rived from the tight binding model ([T). It is desciibed by the 
Dirac theory around the AV, point as 
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with /ill — AsoCTz + aAR2 {kyCTx — k^o-y), where Tq is the 
Pauli matiix of the sublattice pseudospin, tjf = ^at is the 
Fermi velocity, and a = 3.86A is the lattice constant. 
The Hamiltonian H^ explicitly reads 
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in the basis {iJjAt, 'ipBt^ V'a;, V^b;}*, where k± = k^ ± iky, 
and the diagonal elements are 

E{s„t,)=Xsos,t,+eE,t,+Ms„ (4) 

with the spin s^ = ±1 and the sublattice pseudospin t^ = ±1. 
They are not good quantum numbers in general. However, 
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FIG. 4: (Color online) (a) Berry curvature, (b) spin, and (c) band 
structure of a QAH state calculated based on formula (|6) in the Dirac 
theory. Gray curves represent the energy spectrum of a nanoribbon, 
which are identical to Fig|3] Spins rotate by the Rashba interaction 
near the Fermi level, generating a skyrmion spin texture in the mo- 
mentum space. It generates a nontrivial Berry curvature along the 
anticrossing circle whose radius given by (|8). The integration of the 
Berry curvature gives the Chern number C = 2, since there are two 
skyrmions each of which yields C — 1. 



since Ari and Ar2 are very small with respect to the other pa- 
rameters, it is a good approximation to set Ari — Ar2 = in 
most cases. Thus the spin s^ is almost a good quantum num- 
ber in general. An exceptional case occurs when two Dirac 
cones collapse and cross each other, forming a QAH state after 
taking into account the effect of Ar2 7^ 0, as we soon discuss. 

We diagonalize the Hamiltonian ([3]) and obtain four energy 
levels. When two energy levels coincide, the band gap be- 
comes zero, as found in Fig|2l;M,VMP3,VMP2,SVPM). This 
occurs at the K and K' points, where k± — 0. Let us temporar- 
ily neglect Ari because it is very small. Then, the band closes 
when E{sz,tz) = E{s'z,t'z) with dUi. They yield four Unes 
described by E, = ±Xso/i for \M\ < Aso, M = ±Aso for 
\Ez\ < Xso/i, and two Hnes by M = ±{£Xso/Xso)Ez out- 
side the square. They are illustrated by dotted lines in Fig{T] 
These lines are modified by the nonzero effect of Ari, but the 
modification is too small to be recognized in FiglT] See also 
(|5]i for the typical order of correction. 

The Hamiltonian can be diagonalized analytically in some 
cases. First, along the i^^-axis in the phase diagram [FiglT], 
we have already demonstrated!?] that a topological phase 
transition occurs along the iJ^-axis from the QSH insulator 
[FigElQSH)] to the band insulator [FigEtBI)]. The critical 
point is given by 
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where we have set Xri{Ez) = aE^ with a - 
that the effect of Ari is negligible, (ptjO)^ 
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lO-^A. Note 
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SVPM realizes at the critical point, where helical currents 
flow in the bulk. 

Second, along the M-axis, the first Rashba interaction van- 
ishes (Ari = 0), and the energy spectrum reads 



conduction and valence bands with the same spin is given by 



£^±Ja^Xl,k^ 



M 



A|o + fi2^;2fc2 



(6) 



We study a topological phase transition along the Af-axis 
based on this formula (|6]l. When M — 0, there are two 
spin-degenerate Dirac cones for conduction and valence bands 
with a gap between them [Fig|2QSHl)]. As M increases, the 
spin-up (spin-down) Dirac cones are pushed upward (down- 
ward) [FigiaQSH2)]. When \M\ < Aso (l + a^A^j/^^Wp)' 
the band gap is given as A — \M — sAso| at fc = 0, and it 
closes at M = sAso: This is a topological phase transition 
point [Fig|2jM)]. Let us temporally assume Ar2 = 0. Then, 
as I A/ 1 increases further, the two Dirac cones cross each other 
making a circle around each K^ point. Actually, the Rashba 
interaction (Ar2 7^ 0) mixes up and down spins, turning the 
crossing points into the anticrossing points, and opens a gap 
to form the QAH insulating state [Fig|2QAH) and Fig.3]. 
When \M\ > Aso (l + o^Arj/^i^wI), the gap is given by 
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at 



A£± = 2y ft2„2fc2 + ^lE, T Aso)' (10) 

for r]Sz = ±1, and this is independent of M. Thus, the differ- 
ence is smaller for the up-spin Dirac cones at the K point, but 
this is opposite at the K' point. 

We finally determine the phase boundary. It is determined 
as a boundary between insulating and metallic states. The 
Chem and spin-Chem numbers are quantized in insulating 
states, while they are ill-defined in metallic states. As we have 
seen, each Dirac cone moves upward or downward oppositely 
at the K and K' points. Because of this phenomenon the sys- 
tem can become metallic though the gap is open both at the 
K and K' points. This is the VPM state. It occurs when one 
valence band crosses the Fermi level. The condition yields 
four heavy lines {M ~ ±A.so ± iE^) in the phase diagram 
(FiglD- On the other hand, the gap formula dTji determines the 
boundary between the QAH phase and the VPM phase, which 
are the parabolic curves in the phase diagram (FigJT). In pass- 
ing we comment that the VPM phase is metallic in nature and 
does not have mobility gap. Thus the transition from insulator 
to VPM might accompany a mobility gap closing. 
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We present the energy spectrum (|6) and the Berry curva- 
ture calculated by using the corresponding wave function at 
M — 2Aso in FigS As explained there, spins rotates across 
the anticrossing point, generating a skyrmion spin texture in 
the momentum space. This is consistent with the previous 
study for graphene lll4ll . The radius of the anticrossing circles 
is given by (O. We comment that the gap (|7]) is of the order of 
/ieV when AI is of the order of meV. 

We now examine a point in the phase diagram such that 
MEz 7^ 0. In all regions where the effects of Ari and Ar2 are 
negligible, the energy spectrum is derived as 
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The effect of E^ is to change the mass of the Dirac electron. 
Let us increase E^ from E^ = at a fixed value of M. The 
mass decreases (increases) for the Dirac cone characterized by 
775^ = +1 (rjSz = —1) until £'2 — Aso/^, but the behavior be- 
comes opposite after Ez = Aso/^. As a result the tip of each 
Dirac cone is pushed either downward or upward as indicated 
in Fig|2] Consequently the valley symmetry is broken. Note 
that the energy difference at each momentum hk between the 
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